Ac enhancement from strongly coupled quark-gluon plasma 
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We propose the enhancement of Ac as a novel QGP signal in heavy ion collisions at RHIC and LHC. 
Assuming a stable bound diquark state in the strongly coupled QGP near the critical temperature, 
we argue that the direct two-body collision between a c quark and a [ud] diquark would lead to an 
enhanced Ac production in comparison with the normal three-body collision among independent c, u 
and d quarks. In the coalescence model, we find that the Ac/D yield ratio is enhanced substantially 
due to the diquark correlation. 
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The quark-gluon plasma (QGP) is one of the most ac- 
tively pursued subjects in strong interaction physics. Re- 
cent experimental and theoretical studies have revealed 
intriguing properties of QGP. These range from the real- 
ization of the perfect fluid behavior [H and consequently 
of the strong coupling nature of QGP to the lattice 
findings suggesting that the cc bound state could sur- 
vive up to temperatures well above the critical temper- 
ature Tq 0, 0]- Thus a new picture of QGP with non- 
trivial correlations has emerged, and this state is nowa- 
days called the strongly coupled QGP (sQGP) [3,i,il. 

Although the lattice QCD indicates the absence of qs 
bound states in QGP 0, the color non-singlet qq, qg, gg 
bound states may exist if there are attractive channels 
between the constituents Among them, the diquark 
qq with color multiplets and 6c is the simplest one 
[8|. In flavor SU(3)f, the 3c diquarks are classified into 
a scalar with 3f ("good" diquark) and a vector with 6f 
("bad" diquark) [8|, which are in the attractive and re- 
pulsive channels, respectively, according to the one-gluon 
exchange (OGE) picture. Concerning 6c, the diquarks 
belong to either the scalar with 6f or the vector with 
Sf. Although the latter is in an attractive channel, its 
strength is only one-sixth of that in the attractive 3c 
state according to the OGE. Recent lattice calculations 
support the diquark correlation in vacuum . 

The diquark not only has relevance to the color super- 
conductivity at high density and to the Bose-Einstein 
Condensate (BEG) in the strong coupling regime [ll| 



but also is an interesting object for understanding heavy 
baryons Qqq with one heavy quark Q = c, b and two light 
quarks q = u, d [l^ . Generally, in the heavy quark mass 
limit, the li ght quarks arc almost decoupled from the 
heavy quark [l2l | , as the Qq color-spin interaction is sup- 
pressed in OGE by the heavy quark mass |13j and in the 
instanton model by the small coupling between the heavy 
quark and the instanton vacuum as a result of the absence 
of zero-energy modes 14|. Therefore, Qqq baryons only 



in exotic mesons 
In this viewpoint. 



in models that treat the triquark 
as a color non-singlet state 12, 
Ac (Af,) can be regarded as an ideal two-body system 
composed of the c (6) quark and the [ud] diquark, i.e., 
c[ud] {b[ud\). In contrast, the A with an s quark may not 
be such a simple quark-diquark system, because SU(3)f 
symmetry allows interactions among the three quarks. 

In heavy ion collisions, open and hidden charmed 
hadrons are interesting observables for studying the QGP 
21| . particularly at LHC as an appreciable number of cc 



pairs is expected to be produced. Moreover, the ALICE 
detector at LHC is designed to measure charmed parti- 
cles with enhanced vertex tracking system, which has a 
spatial resolution of 12 fj.ni with the best precision [2^ . 
For the planned upgrade of STAR and PHENIX detec- 
tors at RHIC, additional vertex detectors will be added 
to achieve direct vertex reconstruction of charmed par- 
ticles [2^ 24 1 . With such precise measurement of the 



have a strong correlation within the light quark sector as 
in the conventional diquark model 15 , 3, [l3, 0, 13 and 



vertices, even the measurement of open and hidden bot- 
tomed hadrons is possible. 

The existence of diquark correlations in QGP can be 
probed by studying their effects on Ac production in rcl- 
ativistic heavy ion collisions [1^. One of the important 
findings RHIC is that a new hadronization mechanism, 
based on the coalescence of constituent quarks, is oper- 
ative in heavy ion collisions [2^, [l^, [28, 29., Here, 
instead of fragmentation, hadronization takes place by 
the recombination of partons in QGP or by their colli- 
sions into final hadrons. The coalescence model has been 
quite successful in describing the pion and proton trans- 
verse momentum spectra at intermediate momenta as 
well as at low momenta if resonances are included (3l| . It 
also gives a natural account for the observed constituent 
quark number scaling of hadron elliptic flows (3^ and 
the large elliptic flow of charmed mesons [sl]. In such a 
picture, the Ac is formed from the three-body collisions 
among the c, u and d quarks at the critical temperature of 
QGP. If there are strong diquark correlations in QGP at 
this temperature, then Ac could be additionally formed 
from the two-body collisions between the c quark and 



2 



the [ud] diquark. Here, the diquark structure in Ac is es- 
sential to the direct two-body production of Ac, because 
additional process is needed to break up the diquark cor- 
relation if the diquark is absent inside Ac. Since the 
two-body collision generally dominates over multi-body 
collisions we thus expect an enhanced Ac yield in heavy 
ion collisions if there are diquark correlations, and this 
could be a new signal for the search of the QGP. 

The binding energy of the lightest scalar 3t and 3c 
[ud\ diquark can be estimated using a simplified con- 
stituent quark model based on the color-spin interaction. 
This model has been shown to describe very well the 
mass differences between various hadrons including the 
charmed ones [3^ . In vacuum, the color-spin interac- 
tion between two quarks gives the [ud] diquark mass as 
™[ud] = iTT'u + "md ~ Csu • Sd/niumi with the quark mass 
mu = TOd = 0.3 GeV, the spin operator Si {i = u, d), and 
a constant C/ml = 0.193 GeV fitted to the iV - A split- 
ting 3J]. The color-spin interaction effectively contains 
the non-perturbative dynamics in vacuum, and hence 
gives the maximum binding energy 0.145 GeV of the di- 
quark in the strong coupling limit. The diquark mass 
is, however, expected to increase in QGP where the cou- 
pling would become smaller than that in vacuum. In 
the analysis of Ref. 0] , the zero binding of diquark occurs 
slightly above Tc. Since the strength of the color-spin 
interaction is of the same order as the critical tempera- 
ture Tc ~ 0.17 GeV, the diquark correlation could still 
be present near Tc. Therefore, we use the diquark mass 
ranging from TO[ud] = 0.455 GeV for the maximum bind- 
ing to TO[ud] =0.6 GeV for the threshold. 

For the dynamics of heavy ion collisions, we follow 
the expanding fire-cylinder model, which leads to the 
volume Vc — 1000 fm'^ in central Au-|-Au collisions at 
VsiiT = 200 GeV [13 and Vc ^ 2700 fm^ in cen- 
tral Pb-hPb collisions at V^nn = 5.5 TeV At 
Tc = 0.175 GeV, the equilibrium light quark numbers in 
QGP are TVu = iVd ~ 245 ^ and 662 [si] in collisions at 
RHIC and LHC, respectively, all in one unit of midrapid- 
ity. The equilibrium diquark numbers at RHIC and LHC 
for this temperature are estimated as N^nd] — 77 and 208, 
respectively, for mjuj] = 0.455 GeV, and A^[ud] — 44 and 
119, respectively, for mju^] — 0.6 GeV. For the charm 
quark number at the phase transition temperature, we 
take it to be the same as that produced from the initial 
hard scattering of nucleons in the colliding nuclei, and 
their numbers in one unit of midrapidity are Nc ~ 3 and 
20, respectively, at RHIC ^ and LHC We thus 

neglect charm production from QGP, which is unimpor- 
tant at RHIC but could be significant at LHC if the initial 
temperature of QGP is high The charm quarks are 
assumed to reach thermal equilibrium in QGP, and this 
is consistent with the observed large elliptic flow of the 
electrons from the decay of charmed mesons in heavy ion 
collisions at RHIC [s^, H^, which requires that charm 
quarks interact strongly in QGP and are thus likely to 



reach thermal equilibrium [3^, 4^, 41 1 . 

For coalescence of c quarks with independent or uncor- 
related u and d quarks in QGP, the contribution to the 
number of produced Ac is given by 3^, 4^ 



^A°fcud} - 3Ac(cud) 



x/r(2;i--2^n;Pl--P«), 



p, ■ dcr,d^p. 



(1) 



where (7A^(cud) = 2 x 1/3'^ x 1/2'^ = 1/108 is the color- 
spin-isospin factor for the three quarks to form Ac, and 
da denotes an element of a space-like hypersurface of 
QGP at hadronization. Following Ref. [sl], we adopt the 
u and d as well as the c quark momentum distribution 
function fq{x,p) with Bjorken correlation between the 
space-time rapidity and the momentum-energy rapidity, 
and the Ac Wigner distribution function ^cud)(^''P) ~ 

8^ exp(- X^Li yf ~ Y^'i=i^i'^i)' where the relative 
coordinates and momenta are related to the quark 
coordinates and momenta by the Jacobian trans- 
formations defined in Eqs. (7) and (8) of Ref. [s^. The 
width parameter ai in the Wigner function is related to 
the oscillator frequency to hy at = 1/ ^JiMuj with the re- 
duced masses [ii defined in Eq. (9) of Ref. [35|]. Neglecting 
the transverse flow as well as using the non-relativistic 
approximation, we obtain (35l |: 



Arcoal 
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ffAc(cud) 
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2\3/2 



yc(l + 2M«Tca2) 



■(2) 



We note that the Wigner function of Ac used in the above 
does not take into account the [ud] diquark correlation. 
This correlation would reduce the width parameter for 
the relative wave function of u and d quarks. Because the 
number of produced Ac is proportional to the third power 
of the width parameter, treating u and d as independent 
quarks in the Ac thus gives an upper bound for the yield 
of Ac from the coalescence of three independent c, u and 
d quarks in QGP. 

The contribution from coalescence of c quarks with [ud] 
diquarks in QGP to the number of Ac can be obtained 
by setting n = 2 in Eq. ([1]), and replacing the Wigner 
functionof Ac by /^(c[,d])(a^;P) = 8 exp(-yV(T2[^d] " 
'*^'^c[ud])' where y and k are the relative coordinate and 
momentum for the two-body c[u(i] system, and (Tc[ud] = 
1/V^c[ud]^^ with /^c[ud] = TOc"i[ud]/("ic + m,[ud])- Then 
the result is 
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A^c [ud] ) - 5Ae (c [ud] ) A^c A''[ud] , r,„ ■ V ^2 ^ 

l/c(l + 2Aic[ud]JcCrc[ud]i 

with gA^(c[ud]) = 2 X 1/3^ x 1/2 = 1/9. Contrary to the 
coalescence of independent c, u and d quarks from QGP, 
where the \ud\ diquark substructure of Ac is neglected, 
it is here considered as a single entity as assumed for 
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FIG. 1: The yield ratios Ac/D° as functions of temperature. 



the [ud] diquark in QGP. The effect of finite structure 
of the [ud] diquark in both QGP and Ac is expected to 
reduce the yield of Ac in comparison to that obtained 
from Eq. ([3]). The latter thus also gives an upper bound 
for Ac production in the diquark picture. 

The total yield of Ac is given by the sum of above two 
contributions, i.e., N^^;' = iV™(iud) + ^Mcfudl)- The Ac 
yield can be compared with the D meson yield, which is 
not affected by the ud diquark correlation and is deter- 
mined by an equation similar to Eq. ([3]) using instead the 
statistical factor (?do = 1 x 1/3^ x 1/2^ — 1/36, the u 
quark number iVu, the reduced mass fi^u ~ tIcTOu/ {fnc + 
niu), and the oscillator constant CTcu = l/^/icu^^- 

Using the osciUator frequency w = 0.3 GeV, deter- 
mined from the size )ch — 0.124 fm^ of the Df{cs) 
meson based on the light-front quark model [i^, the 



resulting yield ratio Nl°J'^/N^T^ i^c/D°), which is the 
same in heavy ion collisions at RHIC and LHC, is plot- 
ted in Fig. [1] as a function of the hadronization tempera- 
ture. It is seen that the yield ratio Ac/D^ at Tc = 0.175 
GeV is ~ 0.11 without diquarks in QGP (dashed line), 
and it increases to ~ 0.44 in the presence of the diquark 
[ud] with mass mju^] = 0.6 GeV (thin solid line), corre- 
sponding to a loosely bound state which can hardly exist 
near Tq. If the diquark mass has the minimum value 
m[ud] = 0.455 GeV (bold solid line), the Kc/D^ ratio be- 
comes even larger and has a value of about 0.89 at Tc. 
Therefore, the diquarks in QGP raise the Ac/D^ ratio by 
about a factor 4 - 8 in comparison with the case without 
diquarks. For a QGP with a finite baryon chemical po- 
tential, there are more quarks than antiquarks and this 
would lead to a reduction of the yield and an in- 
crease of the Kc/D^ ratio. The expected enhancement of 
Ac/D'^ ratio is still seen if we change the size of Ac by a 
reasonable amount. Because the enhanced Ac/D^ ratio 
is robust with respect to changes in the binding energy of 
the diquark and the size of Ac, more sophisticated models 
are not expected to modify our conclusion. 

^■^From a simple application of the statistical hadroniza- 



tion model 4J], the yield ratio Ac/I?" is roughly esti- 
mated as 2(mA^/m_Do)'^/^ exp(— (mA^ — m£)o)/rc) ~ 0.24 
at the hadronization temperature of Tc = 0.175 GeV 
[s^l- Although this value is a factor of two larger than 
that from the above three-body coalescence of indepen- 
dent c, u, and d quarks, it is smaller than the case that 
includes diquarks. A smaller production ratio is also ob- 
served in elementary processes. In p-l-p collisions, 1630 
Ac's and 10210 D^'s have been measured by SELEX at 
Fermi Lab, and this gives a yield ratio Ac/D^ ~ 0.159 
[4^. In inclusive decay processes of a B meson, the ra- 
tios are is A^/ ~ 0.03 and A^/ D- ~ 0.14 from the 
measured fractions: 79% of D^X and 2.8% of A'X in 
the B+ decay, and 36.9% of D' X and 5.0 % of AjX 
in the B''^ decay, with arbitrary hadrons X [i^. Since 
these experimental ratios include Ac and from decays 
of charmed resonances, a more quantitative comparison 
requires the inclusion of the resonances contribution in 
both the statistical and the coalescence model 



47|. 



The Ac produced in QGP may change into a D me- 
son in the hadronic phase due to collisions such as 
AcTT ND{D*). With the pion threshold momentum 
Pth — 0.43 GeV in the AcTt — > ND process, which is 
larger than the typical energy scale Tc, the conversion 
time due to this process is estimated as l/r = Fth = 
3 crn(p)d'^p/(27r)'^, with 3 being the isospin factor, a 

the cross section, and n{p) ~ exp(— i/p^ -I- m'^/T). With 
(7 = 5 mb as a reasonable value suggested from the J/ip 
dissociation [ist and T = Tq for simplicity, we obtain 
T ~ 17.8 fm, which is comparable^ with the lifetime of 
the hadronic phase tn — 10 fm '4^, leading to a sup- 
pression factor e"*"/"^ ~ 0.57 for the Ac yield. Since the 
temperature in the hadronic phase is lower than Tc , the 
actual suppression factor will be closer to one. Therefore, 
the Ac enhancement is expected to survive the hadronic 
processes. 

In summary, assuming the existence of stable bound 
diquarks in the strongly coupled QGP, we have discussed 
the enhancement of the Ac yield in heavy ion collisions, 
which is induced by the two-body collision between the c 
quark and the [ud] diquark. The Ac enhancement would 
open a new way to find the existence of QGP in heavy ion 
collisions and also provide an experimental tool to probe 
the diquark correlation in QGP. This would, in turn, con- 
firm the diquark structure in heavy baryons with a single 
heavy quark. It is interesting to note that the observed 
suppression of the D meson yield at RHIC 37, 3^ could 
be par tially a consequence of the enhanced production of 

Our study can be straightforwardly extended to At 
production. Using the bottom quark production cross 
sections predicted from the perturbative QCD for p-l-p 
collisions at RHIC [s^ and LHC [lH, we estimate the 
bottom quark numbers in one unit of midrapidity for cor- 
responding heavy ion collisions to be ~ 0.02 and ~ 0.8, 
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respectively. This leads to Ab/B° ratios of 0.098, 0.38 
and 0.82 for the three scenarios of no diquark correlation 
and diquark masses of 0.6 and 0.45 GeV, respectively. 
As in the case of A^, the diquark correlation gives rise 
to a large enhancement in the Ajj/Bq ratio in heavy ion 
collisions. Although the yield of A;, in these collisions is 
much smaller than that of Ac, its enhancement is a bet- 
ter signal to be detected for QGP as the diquark picture 
would be more valid than in Ac, and its much longer life- 
time (t ~ 372/im) than that of Ac (r ~ 62/im) will also 
facilitate its detection. To study the enhancement of Ac 
and At, production is thus an interesting and challenging 
subject at RHIC and LHC. 
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